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A  variational  higher-order  theory  has  been  developed  for  representing  the  bending  and 
stretching  of  linearly  elastic  orthotropic  beams  which  include  the  deformations  due  to  transverse 
shearing  and  stretching  of  the  transverse  normal.  The  theory  assumes  a  linear  distribution  for  the 
longitudinal  displacement  and  a  parabolic  variation  of  the  transverse  displacement  across  the 
thickness.  Independent  expansions  are  also  introduced  in  order  to  represent  the  through-thickness 
displacement  gradients  by  requiring  least-square  compatibility  for  the  transverse  strains  and  the 
exact  stress  boundary  conditions  at  the  top/bottom  beam  surfaces.  The  theory  is  shown  to  be  well 
suited  for  finite  element  development  by  requiring  simple  C°-  and  Cl-continuous  displacement 
interpolation  fields.  Computational  utility  of  the  theory  is  demonstrated  by  formulating  a  simple 
two-node  stretching-bending  finite  element.  Both  analytic  and  finite  element  procedures  are 
applied  to  a  simple  bending  problem  and  compared  to  an  exact  elasticity  solution.  It  is  shown 
that  the  inclusion  of  the  transverse  normal  deformation  in  the  present  theory  provides  an 
improved  displacement,  strain  and  stress  prediction  capability,  particularly  for  the  analysis  of 
thick-section  beams. 
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cross-sectional  area  of  beam 

inplane  rigidities 

width  of  beam’s  cross-section 

the  class  of  continuous  functions  possessing 

discontinuous  derivatives  at  element  nodes 

the  class  of  continuous  functions  that  are 

discontinuous  at  element  nodes 

elastic  stiffness  coefficients 

bending  rigidities 

elastic  moduli 

consistent  load  vector 

transverse  shear  rigidity 

beam  thickness 

cross-sectional  moment  of  inertia  about  y-axis 

element  stiffness  matrix 

beam  span 

force  resultants 

moment  resultants 

applied  transverse  loads 

top  and  bottom  beam  surfaces 

prescribed  end  tractions 

midplane  displacement  along  x-axis 

Cartesian  displacement  components 

transverse  displacements 

axial  and  transverse  coordinates 

transverse  shear  strain 

variational  operator 

strain  and  curvature  components 

bending  cross-sectional  rotation 

dimensionless  coordinates 

Poisson’s  ratios 

stress  components 

beam  element  length 
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INTRODUCTION 


In  the  simplest  cases  of  beam  bending,  vibration  and  stability,  analytic  solutions  can  be 
obtained  using  the  Bemoulli-Euler  (classical)  or  Timoshenko  beam  theories,  with  the  latter 
accounting  for  the  deformation  due  to  transverse  shear.1,2  Analytic  solutions  can  be  either 
difficult  or  impossible  to  obtain  for  the  many  practical  applications  of  beams  as  reinforcing 
members  in  plate  and  shell  structures  and  whenever  nonlinear  deformation/material  behavior  is 
considered.  In  all  of  these  cases,  the  finite  element  method  enables  the  analyst  to  obtain 
sufficiendy  accurate  approximate  solutions. 

The  principal  benefit  of  the  Timoshenko  beam  modeling3'5  is  the  ability  to  properly  account 
for  transverse  shear  deformation,  the  effect  that  can  be  significant  in  deep  beams  and  those  made 
of  laminated  composites  which  are  known  to  exhibit  relatively  low  transverse  shear  stiffness. 
Since  the  Timoshenko  theory  does  not  account  for  the  deformation  in  the  transverse  normal 
direction,  it  precludes  solutions  to  problems  which  lend  themselves  to  the  three-dimensional  state 
of  stress.  An  example  is  the  composite  beam  subject  to  impact  loading6  or  high-frequency 
vibrational  modes  (i.e.,  short-wavelength  loading). 

Recently,  Tessler7'9  has  developed  a  higher-order  plate  theory  for  linearly  elastic  orthotropic 
plates,  incorporating  the  deformation  effects  due  to  the  inplane  stretching,  transverse  shear  and 
transverse  normal  straining.  The  theory  includes  both  linear  inplane  and  quadratic  transverse 
displacement  expansions  and,  due  to  a  special  form  of  a  variational  statement,  is  particularly 
attractive  for  finite  element  development.  An  extension  of  the  above  theory  to  laminated 
composites  and  the  development  of  a  simple  triangular  plate  element  has  been  presented  by 
Tessler  and  Saether.10'11  In  this  paper,  a  higher-order  theory  for  linearly  elastic  orthotropic  beams, 
a  one-dimensional  analogue  of  the  plate  theory,7'9  is  presented.  Also,  a  simple  two-node  beam 
finite  element  is  formulated  and  its  predictive  capability  is  evaluated  with  respect  to  both  analytic 
and  exact  elasticity  solutions. 

In  the  section  on  Higher-Order  Beam  Theory  which  follows,  the  essential  aspects  of  the 
present  beam  theory  are  discussed.  Its  distinguishing  feature  is  that  the  equations  of  equilibrium 
and  appropriate  boundary  conditions  are  derived  from  a  virtual  work  theorem  which  employs 
both  conventional  (Timoshenko-type)  displacement  variables  and  two  higher-order  transverse 
displacement  variables.  The  former  variables  have  the  highest  spatial  derivatives  of  order  one, 
whereas  the  latter  variables  possess  no  derivatives.  The  major  issue  for  finite  element 
development  is  that  (^continuous  finite  element  approximations  being  used  for  the  Timoshenko 
displacements  and  only  (C1)  1  approximations  for  the  higher-order  displacements.  The 
implication  is  that  the  higher-order  variables  can  be  calculated  as  auxiliary  values  and  condensed 
out  statically  at  the  finite  element  level.  This  results  in  the  element  equilibrium  equations  being 
of  the  same  complexity  as  those  corresponding  to  Timoshenko  theory  elements.  The  advantage 


Whereas  C°  continuous  functions  are  continuous  within  the  element  domain  and  at  the  nodes,  C1  functions  arc  only 
continuous  within  the  element  and  can  be  discontinuous  at  the  nodes.  Both  of  these  classes  of  functions  may  yield 
discontinuous  derivatives  at  the  nodes. 
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of  the  present  methodology  is  that  both  transverse  shear  and  transverse  normal  deformations  are 
represented  where  the  resulting  equations  are  the  same  order  as  in  Timoshenko  theory;  however, 
the  range  of  applicability  of  the  present  theory  allows  for  extention  to  thicker  beams. 

In  the  Two-Node  Beam  Element  Section,  a  simple  two-node  beam  element  is  developed  using 
the  present  theory.  The  anisoparametric  interpolations,4,512  originally  derived  to  eliminate  thin- 
regime  shear  locking  in  Timoshenko  beam  elements,  are  used  for  the  element  kinematic  variables. 
Both  analytic  and  finite  element  solutions  for  a  classical  beam  problem  are  presented  in  the 
Discussion  of  Results  Section  and  the  results  are  compared  with  those  obtained  from  Timoshenko 
and  elasticity  theories. 


HIGHER-ORDER  BEAM  THEORY 

In  order  to  clarify  the  development  of  the  theory,  consider  the  bending  of  an  elastic  orthotropic 
beam  having  a  narrow  rectangular  cross-section  of  width  b,  height  2h,  and  spanning  the  length 
L;  the  beam  is  located  in  the  x-z  Cartesian  frame  with  the  x-axis  (xe  [0,  L])  passing  through  the 
midplane  (refer  to  Figure  1). 


WfTfrm 


Figure  1 .  Beam  sign  convention. 


We  expand  the  longitudinal  and  transverse  displacement  components  (ux  and  uz)  with  respect  to 
the  dimensionless  thickness  coordinate  £=z/he  [-1,  1],  where  u*  has  a  special  parabolic  form: 

uz(x,  z)  *  u(x)+h^Q(x),  ut(x,  z)  =  w(x)+^w,(x)+(42-i.)w2(x)  (1) 

where  ^=0  identifies  the  position  of  the  reference  midplane;  note  that  the  expansion  for  u2  is  such 
that  w(x)  represents  a  weighted-average  transverse  displacement  rather  than  the  midplane  deflec¬ 
tion,  i.e. 
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w(x)  =  JL  f\i  (x,z)(  1  -£2)rfz 

4/z 


The  expansion  coefficients  u(x)  and  0(x)  are  also  defined  as  weighted  averages  according  to 


u(x)  =  JL  f *  mx(x,  z)  dz,  0(r)  =  _L.  f*  u/x,  z)  z  dz 


where  u(x)  is  the  midplane  displacement  along  the  x  axis,  whereas  0(x)  is  the  rotation  of  the 
normal  about  the  y  axis.  Assuming  the  beam  is  made  of  a  linearly  elastic  orthotropic  material 
having  its  principal  material  directions  coincident  with  the  Cartesian  coordinates,  a  two- 
dimensional  Hooke’s  law  can  be  written  as13 


^11  ^13  0  £xx 

=  C, 3  C33  0  £„ 
0  0  C«  Y_ 


where  the  elastic  stiffness  coefficients  Ci}  can  be  expressed  in  terms  of  the  engineering  elastic 
constants  as 


C  -  1  C  = _ 

"  l-v„v  ’  33  1- 


V  V 
V31 v  13 


C\ 3  ^33V13'  ^55  ^13’  hV31  ^3V13 


where  £,  and  E3  are  the  longitudinal  and  transverse  normal  elastic  moduli,  v,3  and  v31  are  the 
major  and  minor  Poisson  ratios,  and  G13  is  the  transverse  shear  modulus. 

The  longitudinal  strain-displacement  relation  has  the  usual  form  written  as 


with  the  beam  strains  given  as 


£  =E„+ZK, 

xx  xO  xO 


ho-  M = {“«.*-  ewJt} 


where  a  comma  ()  denotes  partial  differentiation.  The  transverse  normal  and  transverse  shear 
strains  can  be  expressed  according  to  Reference  7  as 
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(7) 


where  the  beam  transverse  strain  quantities  and  associated  through- thickness  shape  functions  are 

{£,».  Ka,  Yj  *{», W/A.  w2(*VA'.  wM.+ew}  _t) 

!*..  ♦,  ♦„}=  {>,,4(4-7?).  >4(3-?).  {(1-?)} 

The  above  expansions  of  the  transverse  strains  constitute  a  major  departure  from  the  conventional 
displacement  formulation.  In  the  present  theory,  the  transverse  strains  satisfy  compatibility  in  the 
weak  variational  sense  by 

5  J*(£„  *  u,fdz  =  0  -  5  £(Y„  "  "  UJ2  dz  =  °  (8) 

where  5  is  the  displacement  variational  operator.  These  transverse  strains  also  ensure  the  exact 
traction  conditions  at  the  top  and  bottom  beam  surfaces  such  that 

xa(x,±h)=0,  oJjc, ±h)  =  0  (9) 

The  beam  equations  of  equilibrium  together  with  the  natural  boundary  conditions  are  obtained 
from  the  virtual  work  principle.  Neglecting  the  body  forces,  the  variational  statement  may  be 


written  as 


j/cxx5exx  +  a«5e„  +  xJ*(JdAdx  -  <7  h)dxdy  -  J  q'(x)but(x,  -h)dxdy 

+  z)  +  Tjut(0,  z)]dA  -  J \tJux(L,  z)  *  Tjut(L,  z)]dA  =  0 

where  S+  and  S'  denote  the  top  and  bottom  beam  surfaces  that  are  free  of  shear  tractions  and 
subject  to  the  transverse  pressure  loads 

xjx,  h)=  0,  on(x,  h)-q\x)  on  S*  ^ 

-h)  =  0,  a t(x, -h)=q'(x)  on  S' 

where  Ti0  and  TiL  (i=x,  z)  are  the  tractions  prescribed  at  the  two  ends  of  the  beam  (x=0,  L),  and 
A  is  the  beam’s  cross-sectional  area  (see  Figure  1). 

Integrating  Equation  10  over  the  cross-sectional  area  results  in  a  one-dimensional  virtual  work 
statement  written  as 

jf  (N£e*o  +  n£Ezo  +  +  +  Q&Jfa  ~  £  +  45vv2)  +  ^5vv,  (12) 

*  Nju( 0)  +  MJQ(0)  *  Qjw( 0)  -  Nju(L)  -  MJQ(L)  -  QJw(L)  =  0 
where  the  stress  resultants  are  defined  as 
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(13) 


(N,  N,)  =  JV„,  aJM,  (M„  Mt)  =  _[(za=,  *,OJdA,  e, - 
(4  is  the  cross-sectional  area ) 

and  the  applied  normal  tractions,  end  forces  and  moments  are  represented  by 

(<?,.  ?2)  =  (<7 + -<7",  q'  +q~),  (q~vqd  =  (qtb,  q2b)  (14) 

-  Jav  rj<M,  (W^  =  JVv  t><m,  (e„  Qj  -  JV„,  T^XM 


where  7^  and  7^  have  the  same  parabolic  distribution  across  the  beam  thickness  as  the  shear 
stress,  tk.7 

In  Equation  12,  expressions  associated  with  the  arbitrary  variations  8w,  and  5w2  must  vanish 
independently,  thus  yielding  the  two  higher-order  transverse  equilibrium  equations2 

NJh -q~2=  0,  MJh 2  -  4^/5  =0  (15) 

Integrating  Equation  12  by  parts  results  in  the  remaining  equilibrium  equations 

Nw«0,  Mxjc -Qx  =  0,  and  Qxjl  +'q‘x  =  0  (16) 

with  the  natural  boundary  conditions: 

At  x=0: 

(17) 

yv  (0)  =  or  8u(0)  =  0,  M  (0)  =  or  86(0)  =  0,  Qx( 0)  =  0xO  or  8w(0)  =  0 


At  x=L: 

Nx(L)  =  yv^  or  8k(L)  =  0,  Mx(L)  =  or  S6(L)  =  0,  Qx(L)  =  0^  or  5w(L)  =  0 


(18) 


Computing  the  stress  resultants  (Equation  13)  by  applying  both  Hooke’s  law  (Equation  4)  and 
the  strain  expressions  (Equations  6  and  7)  results  in  the  beam  constitutive  relations  which  can 
be  expressed  in  matrix  form  as 


N 

X 

*,3 

0 

0 

0 

'  ' 

e* 

N 

2 

^,3 

^33 

0 

0 

0 

M 

X 

>  = 

0 

0 

0,, 

0.3 

0 

4 

M 

2 

0 

0 

0.3 

£>33 

0 

e. 

_0 

0 

0 

0 

G 

Y»o 

(19) 


2 

At  this  stage  there  is  no  integration  by  parts  since  w,(x)  and  w2(x)  have  no  spatial  derivatives  in  the  variational 
statement  (Equation  12). 
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where  the  AtJ,  Dtj  and  G  beam  rigidities  are  given  as 


^11  ~  ^13  ~  ^33  ^33^’ 


Dn  =  C. I,  D.^—IC. 

11  11  Qf/n  >  »3  Qf  J  1 


where 


For  a  rectangular  cross-section 


k2  =  —  and  /  =  f z2dA 
6  * 


A  =  2/zi>  and  /  =  —h2b 
i  3 


(20a) 


(20b) 


Substituting  the  relations  for  N1  and  M,  from  Equation  19  into  Equation  15  yields  the  transverse 
normal  equilibrium  equations  in  terms  of  displacements  which  in  nondimensional  form  appear 


^  <?2 
_  +v.,u  =  — - 

h  13  ■*  2C,, 


hi  '  56C„ 


which  are  readily  solved  for  w,  and  w2.  Introducing  these  variables  into  Equation  19  and  then 
substituting  the  force  and  moment  resultants  into  Equations  15  results  in  the  remaining  equilibri¬ 
um  equations  written  in  terms  of  displacements  as 


(E\AuX  *  (v,,/^  =  0  (a) 

TO  +  |(v13^,  -  k2AGn(Wjt  +  0)  =  0  (b)  (22) 

[k2AGn(w^  +  0)]^  +  ^  =  0  (c) 

As  expected  from  a  linear  (i.e.,  small  displacement)  theory,  the  stretching  (Equation  22a)  and 
bending  (Equations  22b  and  22c)  equilibrium  equations  are  uncoupled.  Observe  that  by  neglecting 
the  Poisson  effect  (i.e.,  by  setting  v13=0),  Equations  22  can  be  reduced  to  those  of  Timoshenko 
theory.  Note,  even  though  a  higher-order  transverse  displacement  expansion  governs  the  deforma¬ 
tion  in  the  present  theory,  the  boundary  conditions  (Equation  18)  are  the  same  as  in  Timoshenko 
theory. 
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In  the  present  theory,  the  solution  procedure  involves  integrating  Equation  22  subject  to  the 
boundary  conditions  given  by  Equations  17  and  18,  then  substituting  u(x)  and  9(x)  into  Equation 
21  to  obtain  Wj(x)  and  w2(x).  Note,  the  stress  components  are  obtained  in  a  consistent  manner 
from  Hooke’s  law  (Equation  4). 

A  TWO-NODE  BEAM  ELEMENT 

To  demonstrate  the  computational  suitability  of  the  present  theory,  a  simple  two-node  beam 
element  is  derived  directly  from  the  displacement  variational  principle  (Equation  12).  The 
interpolation  requirements  for  the  element  displacement  field  stem  from  Equation  12,  where,  as 
in  Timoshenko  theory,  the  weighted-average  displacements  u(x),  w(x)  and  0(x)  have  derivatives 
that  do  not  exceed  order  one;  therefore,  C°  shape  functions  can  be  used.  Equation  12  contains  no 
spatial  derivatives  of  w^x)  and  w2(x),  therefore  the  variables  need  only  be  C1  continuous  (i.e., 
discontinuous  at  the  element  nodes). 

Now  consider  a  two-node  beam  element  of  unit  width  (b=l),  span  I  and  height  2h,  loaded  at 
the  top  surface  by  a  uniform  distributed  pressure  q+(x)=q.  The  displacement  interpolations  can 
be  expressed  directly  in  terms  of  nodal  coordinates  as4,5  which  in  the  present  notation  are 
superscribed  with  an  f,  i.e., 

C°  -Continuous,  Linear  Functions 

u(x)  =  (1  -iDw'o  +  Tiu‘t ,  0(x)  =  (1  -ti)©^  +  r)0*j 

C°  -Continuous,  Quadratic  Functions 

w(x)  =  (1  -ti)w‘0+ti<  +^1(1  -ti)(e‘0-e\) 

Cl  -Continuous,  Uniform  Function 

w,(x)  =  W\,  W2(x)  =  vr'2  (23c) 

where  ti  =  x/C  €  [0,  1],  and  u\,  Q\  and  w\  (i=0,  1)  denote  nodal  degrees-of-freedom  (dof)-  Note 
that  Equations  23  yield  the  beam  strains  that  are  uniform  across  the  element  span 


(23a) 

(23b) 


r-,  =  2K-»'V-2<e,0-ev 

C a ) 

=  2  ®'i  - e*0) 

C b ) 

(24) 

(e,  k*>  -  Vv  wy h) 

h 

(c) 

Substituting  Equations  24  and  19  into  the  variational  principle  (Equation  12)  yields  the  element 
stiffness  equilibrium  equations  in  terms  of  the  six  engineering  dof  («'i(  w'j,  0',  where  i=  0,  1) 
and  the  two  higher-order  displacements,  W *,  and  W*2.  Because  the  latter  variables  are 
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discontinuous  at  the  nodes,  they  are  readily  eliminated  at  the  element  equilibrium  level  through 
static  condensation  resulting  in  the  following 


W\=h 


2  C 


-iVi-O 


33 


w*2  =  h 


56C33  26  1  0 


\lq 


(25) 


Note  that  the  same  result  can  be  obtained  from  the  exact  equilibrium  equations  (Equation  21)  by 
simply  substituting  the  element  displacement  interpolations  (Equation  23). 


Figure  2.  A  two-node  stretching-bending  beam  element. 

The  resulting  two-node  element  is  capable  of  axial  stretching  and  bending  accounting  for  both 
transverse  shear  and  transverse  normal  deformations.  The  element  has  the  simplest  nodal  pattern 
with  six  engineering  dof  (see  Figure  2). 

The  element  stiffness  matrix  K*  and  the  consistent  load  vector  /  corresponding  to  the  vector 
of  nodal  displacement  dof  {u0,  u„  w0,  w„  0O,  0t}*  have  the  form 

Stiffness  Matrix 

ku  kn  0  0  0 

k{2  k22  0  0  0 

0  0  *33  *34  *35 

o  o  kM  ku  k45 

0  0  k35  kA5  k55 

0  0  *36  *46  *56 


0 

0 

*36 

*46 

*56 

*66 


(26) 
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where 


2  hk2Gn 

( 


k  —  k 

j5  *36 


45 


-If  —  ^  U  If  —  lr  —  ^  I 
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Consistent  Load  Vector 


r  -  (/,.  4  4  4  -/3)r 


(26a) 


(27) 


where 

/ia<7V,A  /2=^C-  /,=^(iA2v13“1Lfi2)  (27b) 

The  above  element  stiffness  matrix  (Equation  26)  is  identical  to  that  of  the  two-node  element 
derived  from  Timoshenko  theory;4  however,  the  consistent  load  vector  (Equation  27)  involves  the 
Poisson  ratio  v13  which  does  not  appear  in  the  Timoshenko  element.  (Note  that  in  Reference  4, 
the  axial  deformations  have  not  been  considered.)  Since  v13  appears  in  /„  it  follows  that  the 
element  is  capable  of  predicting  the  deformation  of  the  midplane  x-axis  due  to  the  transverse  load 
q,  which  is  consistent  with  the  theory  of  elasticity  but  is  not  accounted  for  in  Timoshenko  theory. 
As  in  Timoshenko  theory,  the  behavior  of  the  element  in  the  thin  regime  (as  h->0)  is  governed 
by  the  ability  of  the  kinematic  field  to  accommodate  the  Kirchhoff  constraint  of  the  vanishing 
transverse  shear  as  follows 

Y»o  =  w  o> ' -j(®'o +  e\)  -*  0  (28) 


Moreover,  in  this  theory,  the  inextensibility  of  the  transverse  normal  fiber  is  also  enforced  in  the 
thin  limit 

(w,,  w2)  — ►  0  (29) 

In  the  Kirchhoff  constraint  (Equation  28),  the  w\  and  0‘j  (i=  0,  1)  dof  are  consistently 
balanced.  This  is  generally  sufficient  to  guarantee  locking-free  element  performance  in  the  thin 
regime.  However,  shear  locking  can  still  take  place  when  excessive  kinematic  restraints  are 
enforced  on  a  single  element,  such  as  in  the  case  of  a  fixed-pinned  beam;9  e.g.,  when  the 
displacement  boundary  conditions  are  set  w’pO'^w',^,  giving  rise  to  0*, — >0  from  Equation 
28,  which  in  turn  yields  zero  bending  curvature  kx0-»0  according  to  Equation  24,  or  what  is 
known  as  shear  locking.  This  pathological  case,  however,  is  readily  resolved  when  two  or  more 
elements  are  used  in  the  finite  element  discretization.  To  further  improve  element  performance 
and  to  eliminate  shear  locking  from  the  modeling  case  just  described,  a  shear  relaxation 
parameter  (a  device  which  consistently  relaxes  the  Kirchhoff  constraint  at  the  element  level)  can 
be  effectively  employed.5,9 
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DISCUSSION  OF  RESULTS 


To  assess  the  predictive  capability  of  the  higher-order  beam  theory  and  to  examine  the 
behavior  of  the  proposed  two-node  beam  element,  a  simple  bending  problem  involving  a  simply- 
supported,  rectangular  cross-section  isotropic  beam  under  uniform  loading  was  considered.  This 
problem,  having  an  exact  two-dimensional  elasticity  solution,14  encompasses  the  effects  of  both 
transverse  shear  and  transverse  normal  deformations;  therefore,  representing  a  means  for 
evaluating  our  higher-order  beam  theory. 

The  analytic  solution  according  to  the  present  theory  is  obtained  by  introducing  the  uniformly 
distributed  loading  q  (applied  at  2;=1)  into  the  differential  equations  of  equilibrium  (Equation  22), 
which  are  then  solved  for  the  weighted-average  displacement  variables  u(x),  w(x),  and  0(x) 
satisfying  the  boundary  conditions 

m(0)  =  w(0)  =  Mx(0)  =  0 
w(L)=N(L)=M(L)=0 

The  beam  is  taken  to  be  of  unit  width  (b=l),  thickness  2h,  and  it  spans  xs  [0,  L]. 
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Figure  3.  Maximum  midplane  deflection  versus  span-to-thickness  ratio. 

In  Figures  3  and  4,  the  analytic  solutions  of  the  present  higher-order  beam  theory  are  compared 
with  those  of  the  classical  and  Timoshenko  theories,  as  well  as  the  exact  elasticity  solution. 
Figure  3  depicts  the  maximum  midplane  deflection  ut( L/2,0),  normalized  with  respect  to  the  exact 
elasticity  solution,  which  is  plotted  versus  the  span-to-thickness  ratio  (L/2h).  It  is  seen  that  even 
in  the  extreme  thickness  case  of  L/2h=l,  the  present  theory  underestimates  the  maximum 
midplane  deflection  only  by  about  5%.  Although  the  deflection  curve  due  to  Timoshenko  theory 
is  presented  for  comparison,  it  only  represents  a  weighted-average  deflection  and  not  the 
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Figure  4.  Center-span  curvature  versus  span-to-thickness  ratio. 


midplane  deflection,  where  the  latter  quantity  is  unobtainable  from  the  Timoshenko  theory.  Figure 
4  shows  a  normalized  value  (with  respect  to  the  classical  solution)  of  the  curvature  computed  at 
the  center  the  beam  span  (x=L/2)  versus  L/2h.  Both  the  exact  elasticity  and  the  present  higher- 
order  theory  show  identical  solutions  for  the  entire  range  of  L/2h,  whereas  the  classical  and 
Timoshenko  theories  produce  appreciably  erroneous  results  in  the  thick  regime  (for  L/2h<5).  This 
result  provides  evidence  for  the  acceptability  of  the  higher-order  theory. 
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Figure  6.  Distributions  of  transverse  shear  and  transverse  normal  stresses 
across  thickness  for  deep  beam  (L/2h=2). 


Figures  5  and  6  demonstrate  the  through-thickness  distributions  of  the  stress  components 
a„(IV2,z),  a^x, z)  and  t«(0,z)  for  the  case  of  a  very  deep  beam,  L/2h=2.  Whereas  rather  slight 
differences  are  observed  in  the  prediction  for  the  normal  stress  (Figure  5),  the  transverse  stresses 
obtained  from  the  present  theory  agree  with  those  of  the  exact  solutions  (Figure  6). 

The  results  obtained  by  the  finite  element  discretization,  employing  regular  meshes  with  the 
use  of  the  two-node  beam  element,  are  shown  in  Figures  7  through  9,  where  convergence  studies 
for  the  displacement,  strain  and  stress  variables  for  the  deep  beam  case  (L/2h=2)  are  presented. 
The  errors  were  computed  from  comparisons  with  appropriate  exact  elasticity  solutions. 
Throughout,  the  results  are  either  exact  or  converge  rapidly  to  the  exact  solutions  with  mesh 
refinement. 

Finally,  it  should  be  noted  that  the  present  two-node  element,  having  the  same  stiffness  matrix 
as  its  Timoshenko  theory  counterpart,4  does  not  suffer  from  shear  locking  in  the  thin  regime.  The 
element  can  also  be  used  with  the  shear  relaxation  parameter5  to  further  improve  results  for 
coarsely  discretized  models. 


CONCLUDING  REMARKS 

This  report  presented  a  variational  higher-order  theory  for  the  bending  and  stretching  of  elastic 
orthotropic  beams,  including  both  transverse  shear  and  transverse  normal  deformations.  The 
particular  appeal  of  the  theory  is  that  it  provides  a  displacement  variational  framework  for 
developing  effective  and  computationally  efficient  beam  finite  elements  with  the  ability  to  predict 
accurately  the  through- thickness  distributions  of  all  displacement,  strain  and  stress  components. 
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Figure  7.  Convergence  of  maximum  kinematic  variables  for  deep  beam  (L/2h=2). 


Figure  8.  Convergence  of  maximum  strains  for  deep  beam  (L/2h=2). 


A  simple  two-node  element,  derived  from  the  variational  theorem,  demonstrated  improved 
modeling  capabilities  over  a  comparable  Timoshenko  element,  particularly,  in  the  analysis  of 
thick-section  beams. 
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NUMBER  0F  ELEMENTS 
Figure  9.  Convergence  of  maximum  stresses  for  deep  beam  (L/2h=2). 
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